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Abstract 

We study the decay of interface induced structural and electronic relaxation effects in epitaxial 
LaAlOs/SrTiOs heterostructures. The results are based on first-principles band structure calcula- 
tions for a multilayer configuration with an ultrathin LaAlOs layer sandwiched between bulk-like 
SrTiOs layers. We carry out the structure optimization for the heterointerface and investigate the 
electronic states of the conducting interface layer, which is found to extend over two SrTiOs unit 
cells. The decay of atomic displacements is analyzed as a function of the distance to the interface, 
and the resulting exponential law is evaluated quantitatively. 
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Heterostructures based on perovskite transition metal oxides have attracted great atten- 
tion in the last decades , particularly due to the discovery of extraordinary electronic 
and magnetic properties of the internal interfaces. For example, a quasi two-dimensional 
electron gas (2DEG) with an unexpected high charge carrier density develops at the contact 
between the two band insulators LaAlOa and SrTiOa 3|, |J], despite of sizeable band gaps 
of 5.6eV and 3.2 eV . The LaAlOa/SrTiOa heterointerface therefore has been subject 
to intensive research in recent years 3, S, [s]. From the structure point of view it comprises 
(SrO)°, (Ti02)°, (LaO)+, and (AlOa)" layers, where the electron- doped (Ti02)7(LaO)+ 
contact gives rise to the 2DEG. Moreover, the contact between a Mott insulator and a band 
insulator can likewise show conductivity, as realized in the LaTiOa/SrTiOs heterostructure 
[lo| . In all these cases, the phenomenon of conductivity has to be attributed to electronic 
relaxation, because the local charge distribution and its alterations close to the interface are 
the main ingredients. Local deviation from the bulk crystal structure as well as a transfer of 
charge across the interface, induced by different electrochemical potentials in the component 
materials, hence will play a key role for the formation of conduction states HI]. 

Electronic relaxation usually is accompanied by a fundamental lattice relaxation. The 
influence of structure modifications on the LaTiOs/SrTiOa interface has been stressed by 
Hamann et al. [12|, comparing experimental data to results from electronic structure cal- 
culations. In addition to the initial formation of the conduction layer, the mobility of the 
charge carriers likewise depends critically on structural details, such as the incorporation 
of O- vacancies. For the LaAlOs/SrTiOa interface, the mobility increases strongly with the 
vacancy concentration [ij]. Nevertheless, conductivity is maintained in the stoichiometric 
case. Pulsed laser deposition techniques and molecular beam epitaxy nowadays, in princi- 
ple, make it possible to grow layered structures with a precision of a single unit cell, and 
therefore to create atomically sharp interfaces. Local structural and electronic properties 
of such contacts between two perovskite oxide compounds, in general, differ strikingly from 
the bulk materials, similar to the reconstruction and formation of specific electronic states 
at surfaces [14 1. 

Interface atomic structures often are resolved by x-ray diffraction techniques. A contact 
of LaAlOs and SrTiOs layers, for instance, is characterized by displacements of the anions 



and cations in opposite direction, because of Jahn- Teller distorted TiOe octahedra 15|. 
Data of the Ti^+/Ti^+ mixture in the vicinity of the contact suggest, that we actually have 
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FIG. 1: (Color online) Schematic view of the LaAlOa/SrTiOa heterostructure, in a projection along 
the [100] direction and perpendicular to the interface plane. Arrows show the directions of atomic 
shifts due to the interface lattice relaxation. 



no atomically sharp interface, but a metallic Lai-xSr^^TiOs layer is formed [l6|. Moreover, 
conducting tip atomic force microscopy ensures that the metallicity is confined to a rather 
thin region of :^50nm extension 17|. From the theoretical point of view, first principles 
electronic structure calculations using density functional theory (DFT) are a powerful tool 
for studying a structural reconstruction. Results have been reported in the literature for the 
LaAlOs/SrTiOs 



18, Q, 3 and the LaTiOg/SrTiOs interface 3, 2l(. An elongation of the 



TiOe octahedra is found to be a common feature, in accordance with experimental data. 
Electron doping of Ti ions, in combination with the Jahn- Teller effect, leads to the appear- 
ance of metallicity in a band structure calculation applying the local density approximation 
(LDA). Considering the local electron-electron interaction more accurately by means of the 



LDA+U scheme, Pentcheva and Pickett [2^ have succeeded in obtaining a ferromagnetic 
spin order, tracing back to occupied dxy orbitals in a checkerboard of Ti^+ sites. 

While structural distortions have been investigated in detail for different interfaces, an 
analysis of their decay perpendicular to the contact plane is missing so far, although a rather 
general behaviour is to be expected for a large class of materials. The work of Hamann et 
al. [l^ includes data about the atomic displacements but not a systematic characterization. 



Furthermore, t 
is discussed in 



le variation of the electrostatic potential across the LaAlOg/SrTiOg interface 



18(1 . Beyond these data, investigation of structural details as a function of the 



distance to the interface and, in particular, their interrelations to the electronic states is of 
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special interest for a quantitative description of transport processes. In order to solve these 
questions, our present study deals with multilayer structures consisting of a thin LaAlOa 
domain contacted to a bulk-like SrTiOs domain, which allows us to quantify the decay of 
the structural and electronic relaxation effects within the titanate. We carry out a structure 
optimization for the full heterointerface to give a spacially resolved characterization of the 
electronic states. In particular, our data point at an exponential law for the suppression 
of the interface-induced distortion, for which we evaluate the (structural) screening length 
in detail. Comparison to an electrostatic model then enables us to obtain insights into the 
strength of the electronic screening in perovskite-based transition metal oxides. 

Our findin gs r ely on the generalized gradient approximation (GGA), where we use the 



Wien2k code 



23| with a mixed linear augmented-plane-wave (APW) and APW plus local 



orbitals basis set. The package is particularly suitable for describing interfaces [2j, |25|]. In 
all calculations the charge density is represented by ~27,000 plane waves. In addition, the 
k-space grid has 21 points in the irreducible wedge of the Brillouin zone and the Perdew- 
Burke-Ernzernhof parametrization is used. The basis set for the expansion of the wave 
functions contains the valence states La 6s, 6p, 5d, Sr 5s, 5p, Ti 4s, 4p, 3d, Al 3s, 3p, and 
O 3s, 3p, and the semi core states La 5s, 5p, Sr 4s, ip, Ti 4s, 4p, Al 2p, and O 2s. We 
set up a tetragonal supercell by stacking 20 perovskite unit cells: A layer of two LaAlOs 
unit cells in the center of the supercell is sandwiched between bulk-like SrTiOa layers, which 
both consist of nine unit cells. In total, the supercell contains 42 inequivalent atomic sites. 
In the a6-plane, i.e. the interface plane, and in the c-direction the lattice constant is set to 
3.905 A, i.e. the value of bulk SrTiOa. We deal with an n-type interface with (Ti02)°/LaO"'" 
stacking. 

The atomic structure of the LaAlOs/SrTiOs supercell is illustrated schematically in Fig. 
[1] in a projection along the [100] axis, where the interface plane is oriented perpendicular 
to the stack. For convenience, we have numbered the Ti and Sr sites according to their 
distance to the interface. In Fig. [2] we show the partial Ti 3d DOS calculated for the TiO 
and Ti3 atoms of our heterointerface, respectively. These data are compared to the bulk 
Ti 3d DOS. Because of a charge transfer off the LaAlOs layer and, consequently, electron 
doping of the Ti atoms in the vicinity of the interface, the TiO 3d DOS shifts to lower 
energy, with respect to the bulk DOS, and conduction states are formed. As these states 
are characterized by a remarkable (quasi) two-dimensional dispersion, they are attributed 
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FIG. 2: (Color online) Partial Ti 3d DOS (per atom) for the interface atoms (TiO) and the atoms 
in the 4th Ti02 plane off the interface (Ti3), see Fig. [H The corresponding Ti 3d DOS of bulk 
SrTiOa is shown for comparison. 
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FIG. 3: (Color online) Comparison of partial O 2p DOS data (per atom) for the O atoms in the 
first and third SrTiOa unit cell off the interface and bulk SrTiOs. 

to a 2DEG [9], in accordance with experimental findings. However, for a growing distance 
to the interface the doping amplitude declines rapidly, which is visible in Fig. [2] for the Ti3 
site: conduction states have disappeared. We mention that for the neighbouring Ti2 site a 
small but finite DOS remains at the Fermi energy. Moreover, since even the gross shape of 
the Ti3 DOS is similar to the bulk DOS, interface effects have vanished. 

Conduction states therefore are restricted to a narrow area of only two SrTiOs unit cells, 
i.e. three Ti02 planes. However, due to electronic correlations beyond the GGA which tend 
to increase band gaps, this restriction even is only an upper boundary for the extension of 
the electron gas. Our line of reasoning is supported by the O 2p DOS depicted in Fig.O The 
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FIG. 4: (Color online) Setup of the electrostatic model. 

data again refer to the first and third SrTiOs unit cell off the interface and are compared to 
the bulk DOS. While close to the LaAlOa contact the O electronic structure reveals serious 
deviations from the bulk state, only minor differences are left for the O sites in the third 
unit cell. 

The dependence of the structure relaxation on the distance of an atom from the interface 
can be described by a merely electrostatic model. Assuming that all the atomic shifts are 
directed perpendicular to the contact plane, we decompose the system into Ti, Sr, and 
O chains running in this direction. The equilibrium position A = of any atom then is 
determined by the minimum of the electrostatic potential u{A) due to its nearest neighbours, 

AneouiA) = [d + A)"^ + [d - A)-\ 

where d denotes the nearest neighbour distance. In harmonic approximation, the atomic 
force as function of the displacement from equilibrium therefore is 

= — A. 

Let X denote the distance of an atom from the first atom of its chain at the interface, see 
Fig. m which is subject to the initial displacement A{x = 0) = Aq. We obtain 

-1 A(x) = F+(x)-F_(x), 

where F±{x) = d""^ — {d — A{x =l= d))^"^ represents the forces of the two neighbours acting on 
an atom near the interface. A Taylor expansion up to first order gives 

-4A(x) = 2A{x + d)- 2A{x - d) = 4:dA'{x), 
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FIG. 5: (Color online) Calculated displacements in the fully relaxed heterostructure given as 
a function of the distance to the interface Ti02 plane. The numbering of the Ti02 planes is 
illustrated in Fig. [H 
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FIG. 6: (Color online) Logarithmic plot of the data from Fig. [5l A nearly ideal exponential 
dependence on the distance to the interface is obtained for the magnitudes of the Ti and Sr atomic 
displacements. 

and hence A(x) = Aq exp(— rf^^x). Electronic screening is expected to modify the exponen- 
tial decay according to 

A(x) = Aq exp(— a ■ x), 

where we set c? = 1 and introduce an effective structural screening length A > 1. 

The Ti and O atomic displacements found in our heterostructure after the structure 
optimization are shown in Fig. [5] as functions of the distance to the interface, in each case 
given via the number of the Ti02 plane a site belongs to. Moreover, Fig. [5] includes the 
displacements of the interstitial Sr atoms. The amplitudes of all these shifts are largely 



consistent with our model. Specifically, we observe a nearly perfect exponential decay of 
the displacements with structural screening lengths of A = 1.10 for the Ti, A = 1.29 for 
the Sr, and A = 1.56 for the O sites. Corresponding fit curves are shown in Figs. and El 
where the logarithmic representation of the curves in Fig. [H] stresses that deviations from the 
exponential law are negligible for both Ti and Sr. They are slightly larger for the O sites. 
The full deformation pattern induced by the interface lattice relaxation is shown in Fig. 
[U As the Ti and O atoms are shifted in opposite direction, the experimental Jahn- Teller 



distortion of the TiOg octahedra at the interface is reproduced by our data 

The maximal amplitude of the displacements amounts to Aq = 0.023 A for the Ti and 
Aq = O.OI7A for the O atoms, as realized for the interface Ti02 layer. Furthermore, 
extrapolation of the Sr data to this layer leads to a value of Aq = 0.020 A, which seems to 
reflect a similar interface effect on the different atomic species. However, whereas the Ti 
and Sr sites move off the interface plane, the O sites approach it. For this reason, we have 
a relative shift between the Ti/Sr and the O sublattice, which is (exponentially) suppressed 
in the SrTiOs bulk. Consequently, the structural prerequisites for the formation of a 2DEG 
are fulfilled only in the vicinity of the interface plane. In the third Ti02 layer off the contact 
relaxation effects have decayed to less than 4% of their initial amplitude, which explains 
our previous observation that no conduction states are formed here as well as in any layer 
farther away from the contact. Deviations of the structural screening length from the value 
A = 1, as expected in the electrostatic picture, amount to 10% and 29% for the Ti and Sr 
sites, respectively. They can be attributed to the electronic screening, whereas the quicker 
decay of the O displacements is connected to the potential of the O octahedra to evade 
mechanical stress via distortions. 

In conclusion, we have presented first-principles band structure results for a prototypical 
LaAlOa/SrTiOa heterointerface. In general, the electronic properties of heterostructures are 
closely related to the structural distortions affecting the component materials. By means of 
a detailed structure optimization we have investigated the length scala on which such distor- 
tions are suppressed in the bulk materials. We find that the decay of structural relaxation 
effects is well described by an exponential law, in agreement with electrostatic considera- 
tions. However, the structural screening length A is significantly reduced due to electronic 
screening. We expect that these results are fairly independent of the compounds forming the 
interface, as specialties of the chemical bonding are found to play a less important role for 
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the structural screening. In the LaAlOa/SrTiOs heterostructure the formation of a metaUic 
layer is restricted to a narrow area next to the interface, comprising only two perovskite unit 
cells. 
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